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ABSTRACT 


An adaptive control system operating on submarine depth rate 
with constraints on the maximum values of stern plane deflection 
angle and hull pitch angle is formulated using a linear compensation 
technique. The system is adaptive in the sense that it is designed 
to give near uniform dynamic response over ranges of values of various 
hull geometry and speed parameters. The formulated system is simulated 
on a Reeves Electronic Analog Computer, various adjustments made and 
justified, and results shown in plots of stern plane angle, pitch angle 
and depth rate error. A simple system for predicting pull-out time is 
suggested which, when associated with the proposed system, would 


facilitate rapid, well-damped depth changes. 
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Symbols and Nomenclature 


The nomenclature defined in T & R Bulletin No. 1-5 of the Society 
of Naval Architects and Marine Engineers is used herein where applic- 
able. A right-handed cartesian coordinate system is employed with the 
bin at the center of gravity of the hull, the positive x axis dir- 
ected forward, the positive Z axis directed down normal to the base 
line, and angles measured positive in the counter-clockwise direction. 
Pitch angle and depth are measured relative to a similar system fixed 
in space and oriented with the positive Z axis directed vertically 
down. Quantities measured relative to this fixed system are designated 
by the subscript g. All quantities are non-dimensionalized, as indi- 
cated by the prime, by dividing by appropriate powers of the length of 
the hull and linear velocity, the mass and mass density. 


non-dimensionalizing 


symbol fbrmula definition 
工 
I I! y Moment of inertia of 
d Y = Fer hull about Y axis 
k e k Radius of gyration of 
y y = y hull and added mass of 
L hull about Y axis 
k k? k Lateral added mass 
2 2 = 2 Sr 
“m coefficient 
k' Added moment of inertia 
coefficient 
L L' = | Characteristic length 
of hull 
m mt = Mass of hull 





u! 


w! 











1i 


Derivative of moment 
component with respect 
to angular velocity 
component q 


Derivative of moment 
component with respect 
to angular acceleration 
component å 


Derivative of moment 
component relative to 
normal velocity component 
W 


Derivative of moment 
component relative to 
normal acceleration 


Derivative of moment 
component relative to 
stern plane angle 


Angular velocity com- 
ponent relative to y axis 


Angular acceleration 
component relative to 
y axis 


Velocity of origin of 
body axes relative to 
fixed axes in feet per 
second 


Component of u along Z 
axis of body 


Component of acceleration 
of origin of body axes 
along Z axis 


Derivative of normal force 
component with respect to 
angular acceleration 
component 


Derivative of normal force 
component with respect to 
normal velocity component 


ds 


Os 


oL u? 
25 
L 
Ss rad 
6 rad 
1 
U 
Ve? 
subscripts 


Derivative of normal force 
component with respect to 
stern plane angle 


metacentric height 


stern plane angle 


Pitch angle as measured 
between X axes of body 
and fixed coordinate 
systems 


mass density of water 


time in seconds 


desired or commanded quantity 


actual or output quantity 


relative to fixed coordinate system 


iii 











INTRODUCTION 


A. General Aspects of the Depth Control Problem 


Prior to the closing years of World War II, submarine design features 
were based primarily on surface performance. As a direct result of this 
philosophy, submerged speeds, and hence maneuverability, were slow and 
sluggish. The top speed of a typical boat of this period was only nine 
knots at the half-hour battery rate. For this reason the crude two- 
planesman depth control system generally employed since the very earliest 
days of submarine construction was entirely adequate. 

The improved methods of anti-submarine warfare brought about by 
developments in sonar and radar, the use of carrier-bourne aircraft, 
and the construction of fast anti-submarine escort ships made higher 
submerged speeds imperative. Belatedly, submerged performance became 
the primary consideration. By greatly reducing the amount of super- 
structure, and fairing the bridge and periscope sheers, hull resistance 
submerged was reduced by almost a factor of two. This, coupled with a 
doubling of battery capacity, increased submerged speed by a factor of 
two. These higher speeds indicated clearly the need for faster, more 
accurate control of depth. This need is easily recognized when one 
considers that submarines of this period would collapse at depths little 
in excess of one ship length. Even present-day submarines with their 
much higher speeds must confine their motions to a layer of water only 
a few ship-lengths deep. 


The advent of nuclear propulsion plants and the low resistance, 





body of revolution hull form developed in the last decade have greatly 
magnified this depth control problem. Let us consider the situation con- 
fronting the submarine pilot controlling the stern planes when only stern 
planes are used to control pitch angle and depth. This situation is 


“adequately represented by the following block diagram. 





Zgi commanded depth ds stern plane angle 


Zgo actual depth 8 pitch angle 


Fig I 


From this diagram it is evident that there are three integrations be- 
tween the immediate effect of the operator's action and the desired 
result. If the operator is to close the control loop in an optimum 
manner, he must anticipate the results of these integrations - a matter 
difficult in itself, but greatly complicated by the long lag between 0 
and 5, introduced by the lage hydrodynamic mass and small hydrodynamic 
stiffness of the hull. The operator attempts to overcome this diffi- 
culty by introducing a large initial This gives him a more rapid 
indication of the result of his control action, but due to the lag 
mentioned above, he has no timely indication of when to remove or re- 
verse plane angle. This results in overshooting the desired objective 
- the behavior of a poorly damped system. Thus it is evident that the 
operator is required to perform two functions simultaneously in closing 


the control loop: he must speed up the system response, and he must 





provide adequate damping characteristics. This latter skill requires 

a certain amount of prescience highly dependent not only on the skill 
of the operator, but also on his physical condition. In this sense, 
operator reaction is a random variable introduced into the control pro- 
blem which precludes absolutely consistent performance, 

Two basic approaches have been followed to alleviate or avoid these 
difficulties in depth control. The first, and most apparent, is replace- 
ment of the human operator in the control loop by an automatic plane 
angle computer. The second approach, and the one followed in this in- 
vestigation, involves relocation of the human operator outside of the 
closed-loop portion of the system. 

B. Development of Automatic Depth Control Systems 

The German Navy was the first to practice the new philosophy of 
submarine design. Their type XXI boats constructed just before the end 
of World War II in Europe were capable of about twenty knots submerged. 
These ships clearly indicated the desirability of an automatic depth 
control system. An attempt was made to meet this need. The system 
which evolved was really two very similar sub-systems - one for each 
set of planes.” Plane angles were computed as functions of speed, 
depth, and pitch angle. Depth and pitch angle signals were derived 
from pressure gauges and pendulums respectively. Variable gain settings 
in each of the sub-systems permitted ratioing bow and stern plane angles 
as a function of speed and depth-keeping accuracy requirements. The 
controlled quantity was depth, so that the overall system was basically 
a depth regulator rather than an autopilot capable of radical maneuver- 


ing. 





The United States, profiting from German experience, began modifying 
the top-hamper and battery capacity of their 55-475 class fleet boats in 
the late 1940's. Coincident with this Guppy conversion, development of 
a semi-automatic depth-keeping system was initiated. The first such 
system was built by Askania Regulator Company and installed on the 
U.S.S. Irex (SS-482) in 1948. Although very similar to the German 
system in most respects, it controlled not only depth, but also boat 
angle or pitch angle. With this added feature, it provided more man- 
euvering capability than the German system, ? The first of these systems 
were semi-automatic in that the output was an indicated plane angle 
which the planesman was required to match. Subsequent systems were pro- 
vided with an automatic follow-up device making them fully automatic. 

With the advent of Albacore and submerged speeds in excess of 
twenty knots, these early systems became inadequate. Present systems 
are fully automatic, and similar to aircraft autopilots. The Sperry 
System, installed on most new construction, determines stern and sail 
plane angles from signals in depth, depth rate, pitch angle, pitch rate, 
and pitch acceleration to keep the boat on a digitally computed optimum 
trajectory. 

Ce Developments in Manual Control Systems 
The high submerged speeds and limited depth range of modern sub- 


marines make the poor coordination and large time lags inherent in the 


*The Guppy (Greater Underwater Propulsive Power) conversion involved 
streamlining the superstructure and doubling the battery capacity 

of the SS-475 class boats. These two modifications resulted in a 
200% increase in submerged speed. The resulting ships bear an un- 
canny resemblance to the type XXI boats. 





old two-operator, peso control system intolerable. The first step 
away from the old system was incorporated in Albacore with the substi- 
tution of an aircraft-type control column for the old wheel type de- 
en Two sticks were provided. The main unit is a conventional 
aircraft steering column in which wheel rotation controls rudder posi- 
tion, and fore and aft motion of the entire column controls bow plane 
angle. The second stick controls stern plane angle. In the principal 
mode of operation, stern plane angle command is derived from the bow 
plane signal via a variable gain amplifier, thus providing single- 
operator ratioed control. This is the mode of operation greatly pre- 
ferred by the operators. 

This basic system has proven satisfactory in most respects. Sub- 
sequent developments have been primarily concerned with providing the 
operator with the required feed-back information ~ depth error and boat 
angle - in the most efficient manner. The first major improvement in 
this line was the Combined Instrument Panel which displays all the 
information required by the operator in a cluster of dials arranged 


Sm: 


in an optimum manner from the human engineering point of view. 
second, and rather novel, development is the Contact Analog Display 
developed by the Electric Boat Division of General E In es- 
sence, this is a television screen picture of the path the operator is 
required to follow. A three dimensional effect is provided by a fic- 
ticious horizon and narrowing-with-distance perspective of the pic- 
tured road. This is the only indicator the pilot must concentrate his 
attention upon. His task is simply to keep the white line of the road 


properly positioned on the screen. 


In anticipation of the performance capabilities to be required of 
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future high speed submarines, the Bureau of Ships has initiated an ex- 
tensive study of the over-all submarine control problem. The title of 
this study is the Submarine Integrated Control Program (SUBIC). The 
Contact Analog mentioned above originated in this program. As a part 
of this investigation, Cornell Aeronautical Laboratories conducted a 
study of the data flow and data processing requirements in effective 
submarine soon! One of the major results of this study was an 
indication of the desirability of locating the human operator external 
to a multiple closed-loop system in which the outer loop closure is on 
depth rate. The advantages of this scheme are readily apparent. In 
addition to minimizing the effects of erratic human behavior and sim- 
plifying the operators' problem by avoiding two of the integrations 
mentioned in (A) above, the system is easily made fully automatic by 
substituting a digital computer or programmer in place of the human 
operator. This system is represented by the following block diagram, 


The K's indicated are in general non-linear functions of ship speed and 


2 
Computer 


Zoi 





Fig. II 


the geometrical parameters of the hull. 





DE The System to be Investigated 

The major difficulty in submarine control system design is the lack 
of precise information concerning the values of the various hydroframe 
parameters. At the present time, there is no satisfactory method of es- 
timating these parameter values within acceptable accuracy from the 
geometry of the hull as delineated on the lines drawing. Recourse to 
model and full-scale tests is unavoidable, 1? Further, in order to 
obtain acceptable performance over the speed range of a given hull, the 
characteristics of a control system such as that shown in Fig. II must 
be varied in some pre-determined, generally non-linear manner with 
speed, It is with this problem that this thesis is concerned. 


The basic system under consideration is shown in Fig. III. It is 


evident from this diagram that no provision is made for a steady state 
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Mage LIT 
plane angle in a steady dive. This feature would require an auxiliary 
forward loop from the input to the plane angle summer. This matter is 
ignored in the following investigation because the effect of such a for- 
ward loop is small and decreases with speed roughly as the reciprocal of 
the speed cubed. The effects of this simplification will be a slight 


increase in steady state depth rate error at slow speeds. Also, only 





stern plane control is considered. 
The objective of this investigation is to formulate suitable com- 


pensation functions KÒ E: and K „G that will yield a specified per- 


l cl 02 62 

formance of the closed loop system over a specified speed range of a 
given hull configuration, and at the same time minimize the effects of 
deviations of certain important hull parameters from their specified 
values. In this sense, the proposed system is an adaptive system. 

The problem will be treated on a non-dimensional basis using the 
linearized force and moment equations of a submerged body in the ver- 
tical plane, and standard nomenclature as specified in reference [1]. 
The system will be simulated on a RBAC analog computer to determine its 


transient and steady state response to a step input of depth rate under 


various conditions. 
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PROCEDURE 


A. Scope 


This investigation is concerned with the formulation of a suitable 
closed-loop system for controlling the depth rate of a deeply submerged 
submarine subject to constraints on the maximum values of hull pitch 
angle and plane deflection angle. The principal objective is a system 
whose dynamic response to an ordered change of depth rate is relatively 
independent of variations in ship speed and hull geometry. A linear 
design technique will be employed, and the performance of the resulting 
system analyzed on an electronic analog computer, 

The basic idea underlying the approach to be followed is that 
hydroframe parameters may be considered to fluctuate from hull to hull 
as well as for a given hull about certain normal values. Similarly, 
ship speed hey be considered as a deviation from some normal speed. 

The effects of these deviations is a change in the resonant frequency 

of the hull and hence of the entire control system. If when its para- 
meters have their so-called normal values, the overall system is so 
designed that its resonant frequency is much higher than any frequency 
to which it is normally required to respond, then the change in natural 
frequency brought about by parameter variations will have a small effect 
on system performance. 

The design specifications for the proposed system are as follows: 

l. a non-dimensional natural frequency of the "normal" system 


in the vicinity of 5 with a damping ratio of about .75. 





2. no more than 15% overshoot in pitch angle.* 
3. no more than 10% overshoot in depth rate. 
4. a rise time to 10 pitch angle for large changes of depth 
rate close to 1 second in non-dimensional time. 
The non-dimensional form of the equations of motion of the hull 
are employed since this facilitates separation of the dominant para- 


meters from those of lesser importance. 


B. Equations of Motion of the Hul110 


Omitting the cross-acceleration terms M. w' and SI which are 
usually quite small,” the linearized equations of motion of a submerged 
body in the vertical plane may be re 


(1) (m'-2,') w - AN = (m! + 20) q! = 25 5, 


S 


_V" 1 sr t Ve ' won d 0 도 0 
(2) -Mi w + (1, M) å Mr a'-a' 2; 0 KS 


The coefficients of the acceleration terms w' and q' are inter- 
related through the following equations: 
2 
ca m' 
e Y ) 
(gege t 
Ze kom 


해. 1 -k' Tt 
q Es y 


k' is simply the non-dimensional radius of gyration of the hull about 
the y axis. k, and k' represent the added mass of the hull. They are 
usually taken as the theoretical values applying to the ellipsoid of 
revolution approximated by the hull SE 

*It should be noted that the desired glide angle is implicitly included 
in the commanded depth rate since Z =w- uO. Assuming a steady glide 


with zero normal force, W = o and ZO = - 0 0, In non-dimensional nota- 
tion, Z} =- © S 


10 





Landweber and Johnson of the David Taylor Model Basin have pro- 
posed the following relations among the coefficients of the velocity 
terms for the case under consideration where the only appendage is a 


tail e 


os 
H 


COL 

These relations, although only approximate, give predictions of the 
indicated parameters which are reasonably close to measured values. 

x, may be interpreted as the non-dimensional effective moment arm at 
which the lift due to angle of attack acts. The m' in this relation 

is taken as an approximation to the Munk moment which appears in the 
potential flow about the body. X, may be regarded as the non-dimen- 
sional effective radius through which the angular velocity q is related 
to an equivalent normal velocity, i.e. w' = E oi, (x3) may be inter- 
preted as the effective product of the radius relating angular and 
linear velocities and the moment arm relating moment and force due to 
angle of attack. Although such a product is not necessarily a square, 
it seems nearly so and is so defined for convenienge 주 


The coefficient of 6 in equation (2) may be reduced to a speed 


parameter as follows: a ou 
o 








o Co -- EU bg Eai É 
m' | 5 y? ye 
ark 


This speed parameter Ù is then dependent not only on the ship speed U, 


but also on the metacentric height Ze which is a measure of the ship's 


11 





longitudinal static stability. Incorporation of the minus sign in the 


definition of 1/0 is done since Z, is always a negative quantity for 


B 
ships with positive static stability. 
The formulation is further simplified with the introduction of the 


two quantities 





Z' 
En è > 
3 Zi 
W 
M’ 
A! = 一 3, 
21 
S 


Ki may be regarded as a measure of the control effectiveness or as the 
ratio of the lift force generated by the stern planes to the lift force 
generated by the fixed surfaces of the hull due to angle of attack. X, 
is the non-dimensional effective radius at which the normal force gen- 
erated by the stern plane acts. The negative sign arises from the 
coordinate system employed. 

Making all of the above substitutions in equations (1) and (2) 
after dividing through by m!’ yields the equations of motion of the 


hull with which we shall be concerned. 


(3) CRE NT e E 
= m' m' g 
1 0 ? me 9 nt - IVER. d wm =X 171 
(4) SOI ee Ap JW + (e DO (1 + kt) à za re IR KZ, e 
mi m' U mu g 


AS 
The utility of this formulation lies in the fact that U is the 
only speed-dependent parameter of importance, and the number of geome- 
try-dependent parameters has been greatly reduced. If it is assumed 


that an varies as the first power of the speed, which is not unreason- 


de 





able, then it can be seen that 2 '/m' is independent of speed from the 





following: CN 
Kol“ U 
Zu Mi 
Z = = W 
~ 3 Um 
FAL 
For hulls of normal design, the ranges of variation of Kos k', and 


Ey are very limited. Similarly, the effective moment arms and radii 


1 


values that they may be considered to be of secondary importance in 


… x,» and (a, 1)" are confined to sufficiently limited ranges of 


this investigation. Accordingly, these parameters are assigned the 


following normal values and henceforth treated as constants. 


x! = -0.5 000, TOAS 

zo! = 0.5 KÉ wa 0.933 

Gr EE k' = 0,805 
x! = -0.5 


Then the dynamic characteristics of the hulls to be considered 
AS 
depend upon the three primary parameters U, Zw'/m', and Kot. As 


noted above, the metacentric height Z also enters implicitly as a 


B 


forth major parameter. These primary parameters will be considered to 


fluctuate about assigned normal values as follows: 


AS 


i+ 


E ou 05.02 D.o 
E =.-1.5 > 0.5 
mi 


: = e e 
Kt 0.3 + 0.1 


It is considered that the indicated ranges of Zw' /m' and K 


include most body of revolution hulls existing cy to be constructed. 


23 





AS 
The actual speed range indicated by the range of U depends, of course, 


on the value of metacentric height selected. This figure normally lies 

in the range -O.5 to -2.0 feet. For Ze = -0,5 feet, the indicated range 
nd 

of U corresponds to a speed range of 2.3 to 21 knots. For Zp = —2.0 

feet, this range corresponds to speeds from 4.7 to 41 knots. Again, it 


is considered that the range of U considered includes all combina- 


tions of Ze and U of practical significance. 


C. Formulation of the Proposed System 


By applying Laplace transform operations, in which the operator "s" 
implies differentiation with respect to non-dimensional time , to 
equations (3) and (4), one obtains the following transfer functions for 
the hull. 


(5) els). = (E, s + F,) s 
SA 3 


+ B gé +Cs+D 











As 
3 2 
1 = = 
WE) 2" (a) . Mes N, s Na 
9 (8) As? 十 Be" + Cs + D 
where 
2 
= 1 t 
A (o), VC (1+x') 
1 
E (ki e (l+k') + (1+k,) (x «yé 
RRE A SA 
2 We t cm 
Cn (2) ay)? GES co! | NE ED EE + 1 
m' m' U 
D= -2']1 
Zar men 
m' ye 
= SÉ 0 9 0 
E, = en SI EI A 
m 
o Q t ' 
F= 2' K'| 2, (x, x] Jun ds 
m' m' 
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: 2 
ATA (1) (1+k' ) 





m! 
N za ' t ' 
2 "e Sp ~ (x,' e E = Ka. 
m' m' S 2 8 
N = 1 1 e = 
3 ko an | (x,' x DEL l + i 
1 1 
m m ue 


Equations (5) and (6) are rewritten in standard from as follows: 


SC) e S = K, (s+H)s 





o ls 
S 3 2 
SA dp + J,8 + u 
8) Zt 2 
(8) Zo (s) n K, (8 + E,s+ n 
5 (s) 3 2 
3 s + JS + J, s+ d. 
where 
Ki = E, K, = N 
A A 
N = P. E, = -N, 
E, N 
F 
= -N 
ee 
A 
J, = B d. = D 
À A 
J = C 
L A 


Substituting the assigned values of the secondary parameters 
given above, these coefficients are expressed as functions of the 


primary parameters as follows: 


d ' a ' d 
Kl = 6.16 K, ~ E, = 0.517 Ko Ze 





m' mi 
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"= 1.034 Ep 75-145 


don -3.6 0 ' He (ll) 
zn 2 
m u 
J= 63( 1.93: A a) 1, = 637 A 
U m' we À 
U m 


The following points are now evident concerning the transfer 
functions defined by equations (7) ana (8): 
1. both gains are always negative since Zw' is a negative 
1 

number and K is always positive. i 

2. both are stable transfer functions in that neither has 
a right half plane pole, 

3. (8) is non-minimum phase since E, is negative. 


Using these relations, the system to be investigated is re- 


presented by the block diagram of Fig. IV. 





Fig. IV 


The open-loop relation for the unity feedback system is: 
i K_G 
ON mm (O) vw 19 (5) 
ee cl cl (3) IKE o0 T) KG, (s) 
e (s) 1 02 1 02 
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where Kel and K „G are the compensation functions to be formu- 


Ge Ce 
lated, and 
d (8) = K, (s+N) s 
3 2 
SA Jo Ss + Y, S+ d. 
G, (s) = K (a? + E, s + F,) 
2 2 2 2 


K, (8 + N) 8 


es) = Zi (s) - 2! (s) 


gi go 
Substituting the expressions for K} G} (s) and K, 6, (s) into 
(9), one_ghtains 
1 
(10) D LY K KG (s) (ES + E,s + F 0) 
e (8) = ee re E (s)(s+N)s 
2 1 O Te c2 


Now we choose G , (s) to be of the form 
(8) = e a: 
s(s+N 


in which we note that N is a constant and therefore independent of 
parameter variations, 


With G > (s) of this form, (10) becomes 


(11) 78 (99 si O) (s% E, 


3 2 
8 + ~ KE Je + (J4- KK als +(J -K 


8 + $.) 





e (s 
120) 


and the utility of this selection of G(s) is immediately obvious. By 
suitably choosing Ko) a, and b, the poles of the open-loop relation 
may be conveniently located in the s plane. (As will be indicated 


later, the situation is not quite so straightforward due to saturation 


ind.) 
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Ås stated above, the system will be designed with the primary 
parameters assigned their normal or mid-range values. Inserting these 
numbers in (11) yields 
(12) 2 (a) -.233K 1G 1(s)(s“- 5.7458 - 11.8) 

e (s) 5 +(5.4 + 2.77K ,) só + (3.68 + 2.77K ja) 8 + 
(.382 + 2.77 K_.b) 

For the present, take G (8) equal to unity. Then the zeros 
of the open-loop relations are located at + 7.35 and - 1.61. Since 
Kel must be positive, we are concerned with variations of the negative 
gain factor -.233 Kal as far as stability and performance of the closed- 
loop system are concerned. The locations of the poles of (12) are 
determined by experimenting with the root locus plot of this relation 
to achieve the desired resonant frequency and damping ratio of e 
dominant pole-pair of the closed-loop function. The resulting config- 
uration is shown in Fig. V. It is to be noted that the most negative 
pole of (12) is one of the dominant poles. 

With the poles of (12) located at -1, -4, and -9 as shown, the 
denominator of (12) becomes 

GE eet 49s + 36 = 874 (5.4 + 2.71 또 Ee + (B66 4 2.77 K als + 


(.382 + 2a ob) 


Equating coefficients of the various powers of s yields 


a = 5.26 


b = 4.14 


The phase-magnitude plot of the open-loop relation incorporating 
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these values is shown in Fig. VI as curve A. From this plot it is 


evident that D. should be a lag network in order to increase the low 


1 
frequency gain of the system. Taking G_, (s) to be given by 


G (s) = § +2 
다 Ss + sé 


curve A is modified as shown in curve B of Fig. VI. Applying curve 
Bto a Nichol's chart indicates a required gain factor of 14 db for 
an SÉ of 2.3 db. That is 


20 log .233 K j = 14 


' = 
KÒ 21.4 
The frequency response of the closed-loop system is shown in Fig. VII. 
A steady state error of -.4 db, or 4.59, is indicated in depth rate, 
but this is considered acceptable. 


With Ki H (s) and Ko G s) thus determined, the overall 


ce ( 
system was simulated on a Reeves Electronis Analog Computer. As will 
be indicated, it was necessary to modify G , (s) in order to realize 
a satisfactory step response. The computer set-up is described in 
Appendix A, and the plotted responses in 6 EL 8, and e for a step 


input of EC are shown in Appendix B. 
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III 


RESULTS 


A. General 

The proposed system was simulated on a Reeves Electronic Analog 
Computer as described in Appendix A. A constraint of plus and minus 
25 degrees or .436 radians was placed on the maximum stern plane angle 
permissible. The input ¢o the system was chosen to be a step function. 

Zt = .5 u_l(z) 

This choice was made on the basis that in a steady glide w' should be 
desirably zero making 


2' = -6 
q 


where 6 is the maximum permissible steady pitch angle. This value of 
8 was taken to be -.5 radian or -28.6°. 30° is usually taken as the 
maximum allowable value of ® in medium speed submarines. The actual 
peak value of O indicated in the results of Appendix B are only ap- 
proximate due to the assumption that 8 = sin O. The error introduced 
by this approximation is justified in he sense that using a high 
value of © permitted observing the system's general performance in the 
most extreme condition normally encountered. 
B. Evaluation of Results 

With the system as designed in Chapter II set up on the computer, 
and the primary parameters assigned their normal values, it was observed 
that the overshoot in O and e were of the order of thirty percent. This 
was due to excesSive lag introduced by Gi and saturation in 6 5: Ac- 


cordingly, the locations of the pole-zero pair of the lag network were 
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shifted to -0.1 and -1.0 respectively by changing the setting of po- 
tentiometer number 15. This corresponded to changing the value of 
k in the lag function representation given in Appendix A from l to .5. 
The effect of this adjustment on the phase-magnitude plot of the open- 
loop function and on the closed-loop response are shown in Curve C of 
Figs. VI and VII. After this adjustment, the system performed 'satis- 
factorily. Plotted results are shown in Appendix B., 

The following observations are to be noted concerning these 
results: 


l. the system specifications are met for values of Ze? of 
m! 
-1.5 and -2.0 for all values of U greater than 2.0 and 


all values of Kè considered. 


2. the system does not function satisfactorily for Z z=],0. 


m' 
~ 
At values of U of 1,0 and 2.0, the system either responds very 
slowly and with large steady-state error, or shows very little res- 
ponse at all, especially for low values of Kie This is due directly 
to the small lift force generated by the stern planes at slow speed 
and cannot be taken as a failure of the plane angle control system. 
At these speeds, 8. takes on a maximum value and retains this value 
as would be expected. For this reason values of less than 3.0 are 


not considered in the results, 


el 
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IV 
DISCUSSION OF RESULTS 


1t was noted that the system did not function satisfactorily 
for Z'/m' equal to -1.0. The reason for this becomes apparent 
when the effects of this parameter and of saturation is O, on the 
locations of the open-loop poles of the system are considered. 
Limiting the maximun value of Se is equivalent to limiting the ef- 
It is recalled 


fectiveness of the compensation gains E, and KÒ 


1 2" 


that the open-loop poles of the system were moved further out into 


the left half of the s plane by a suitable choice of Ko The effect 


2º 
of saturation ind, is,therefore, to constrain these poles to lie to 
the right of their desired locations. The effect is aggravated by 


AA 


the tendency of J. as well as Je to approach zero as U increases when 


1 


Z'/m' equals -1.0. For values of Z /m' between zero and -1.0, d. 
approaches a negative number. In the extreme case, two of the poles 
fall between the origin and the zero of the lag network, and these 
become the dominant pole pair. The net result"is a lowering of the 
closed-loop natural frequency and damping for a given gain. Two 
possible remedies present themselves. First, the pre-saturation gain 
Ki may be decreased to keep the system in the linear region of 94 
longer, with a resultant decrease in steady-state accuracy. Fig. XXVI 
of Appendix B shows this result for Luz 9.25 instead of 21.4. A 
second approach lies in increasing the post-saturation gain of the 


hull via Kou This is equivalent to supplying more control surface 


area. Figs. VIII through XIII of Appendix B indicate this tendency. 
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A more detailed discussion of the effects of saturation indicated 
would require a non-linear analysis of the system and will not be 
attempted here. 

The question immediately arises as to why these effects were 
not predominant at the higher values of Ze investigated. The 
necessity of shifting the zero of the lag network to the right men- 
tioned above indicates that these effects were present when Beim: 
equalled -1.5, but a suitable remedy was quite simply achieved. The 


effect is not so pronounced simply because J, and d assume larger 


1 
values and the hull gain is larger for these values of 2 '/m'. 

This failure of the system at low values of Ze is not considered 
a serious limitation on the basic approach presented. For normal hull 


15 


forms, the magnitude of Zn /m' is rarely less than 1.2. The re- 
sults simply confirm the known fact that hulls with little directional 
stability must be provided with larger control surfaces. 

The remainder of the plotted results of Appendix B for values 
of Zt /m' of -1.5 and -2.0 are considered to justify the method pre- 
sented. The behavior of pitch angle, a critical quantity, in the dive 
entry maneuver simulated is considered to be quite satisfactory. The 
error in steady state depth rate is maintained less than eight percent 
in all instances, and the overshoot is negligible. It will be re- 
called that the calculated steady state error of the system as com- 
pensated was four and one-half percent. From Fig. XVII, the steady 
state error of the system as designed is seen to be four percent. Since 


PA 


U appears in the denominator of the expressions for J. and Jos decreas- 


1 
~ 
ing U will move the poles of the open-loop relations to the left and 


es 





thus increase the steady state error, an effect clearly shown in the 
plotted results. e e is true when y is increased. It is 
seen from the plots that this effect is made quite small by the 
compensation technique employed. 

1t will be recalled that the requirement for a steady-state 
plane angle to counteract the metacentric righting moment was neg- 
lected in the system formulation. From the results, it is seen that 
this requirement is almost completely satisfied by the system due to 
the effect of metacentric righting moment on depth rate. For low 
values of Ù, an appreciable steady state plane angle exists but 
there is still a small pitch rate. 

It is noted that the effect on steady state error due to changes 
in hull gain through variations in Ge and 2, /m are smaller than the 
effect of metacentric moment noted above. However, several runs were 


made with Ka adjusted to correct for this change in hull gain. The 


1 
runs were made only for the case U = 5.0, but the results shown in 
Pigs. XXVII and XXVIII indicate that a slight improvement in steady 
state error results. 

In the basic system presented, it was assumed that an infinite 
plane rate, $, was available. This is not the case in practice. The 
compensation method presented, however, can easily be extended to 
handle the case when plane actuator dynamics are included in the KG, 
function. For example, if these dynamics are representable by a 
transfer function of the form nn addition of the term (p s+1) 


p s+1 


in the numerator of both Gi and Go would cancel this function and 


require minor changes in the constants associated with K 1761 and 
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KG ye Again, complete cancellation would be achieved only in the 
linear region of stern plane operation. 

It is becoming common practice to evaluate performance of a 
proposed submarine control system in terms of the comparison between 
the trajectories achieved for step inputs to the system and experi- 
mentally determined minimum time or limit trajectories. Due to time 
and computer capacity limitations, such a direct comparison was not 
carried out in this investigation. lt is maintained, however, that 
on the basis of the indicated behavior of stern plane angle the 
plotted trajectories of depth rate error and pitch angle for the 
cases of interest closely approach those observed in the limit dive 
entry. In the results of Appendix B for the cases in which 20 /m' 
assumes the values of 1.5 and -2.0, 2。 shows at most one complete 
reversal from its initial maximum value before it quickly settles 
out to its steady state value. This behavior is necessarily charac- 


teristic of a minimum time trajectory. More will be said on this 


point in the following chapter. 
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A PROPOSED PULL-OUT TIME PREDICTOR 


The proposed system functions satisfactorily in carrying the ship 
through a dive or rise entry phase and maintaining the required dive 
or rise, but it cannot inform a human operator when to pull out to 
achieve a desired depth. As was indicated, a computer or programmer 

could be substituted for the human operator, but a simpler solution 
would be desirable. In the course of the investigation reported 
herein, a simple analog method of computing pull-out time presented 
itself. The basic idea stems from an optimum controller concept 
proposed by Lieutenant M. C. Davis USN in his doctoral thesis on multi- 
dimensional random processes... This technique is considered worthy 

of note and further investigation as an auxiliary component of the 
proposed system. 

If the optimum trajectory in depth is specified to be one executed 
in minimum time with no overshoot, then the optimum controller applies 
maximum control effort throughout the entry and pull-out Haden Re 
Consider a ship controlled by the proposed system to be in some steady 
dive. In pull-out, the operator could simply return his control col- 
umn to the zero position, commanding zero depth rate, or he could com- 
mand varying amounts of negative depth rate to speed up the pull-out. 
But we have seen that the system presented in this paper approaches 
the performance of the optimum controller itself when the input is the 
desired step change in depth rate. Therefore, commanding a negative 


depth rate in pull-out cannot give a faster response than simpiy 
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returning the control column to zero. The problem, then, is deter- 
mining when the operator should Zero his control column to level out 
at a desired depth. 

Th problem and an indication of a possible solution are shown in 


the following figures. In (a), the ship is shown in some steady dive， 
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as indicated by the solid line, which will carry it through the de- 


sirsd depth Z At successive instants, to the present state of 


gi" 
the system.as determined by speed, depth, depth rate, pitch, and 
pitch rate are introduced as initial conditions into a fast time 
analog of the system, shown in (b). The input to the analog is a 


step function of equal magnitude but opposite sign to the existing 


commanded depth rate, a step function returning commanded depth 


eT 


rate to zero. The computing periods, in real time, are short relative 
to the response time of the hull. The output of the analog is inte- 
grated to yield, in fast time, the trajectory in depth which would 
result if the control column were zeroed at the associated sampling 
instant tj 。 This projected depth response is subtracted from the 
desired depth, represented by a constant quantity, and the difference 
is supplied to a zero-crossing detector. At the first zero-crossing 
of this difference, the indicator is actuated. Now it is noted from 
the dashed representations of the projected trajectories in (a) that 
at the first zero in this difference both the depth error and its 
first derivative are zero. That this will be the case is indicated 
by the plots of depth rate error response in Appendix B. This simp- 
lifies the operator's task since if he zeroes his control column at 
the first indication of a zero crossing he will come out on depth 
with depth rate passing through zero. Then he need only correct for 
the error in depth due to the small overshoot of the depth rate con- 
trol system. Note that pullout time is based on the projected depth 
at peak overshoot in depth rate so that in a dive this error will 


always be on the safe side. 
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VI 


CONCLUSIONS AND RECOMMENDATIONS 


On the basis of the results shown in Appendix B, it is concluded 
that the system presented is entirely feasible and worthy of further 
investigation and implementation. It is noted that this system pos- 
sesses the following very desirable characteristics from the operator's 
viewpoint: 

1. the system adjusts itself to any steady dive or rise 
commanded by the operator in a minimum time. 

2. response becomes faster as speed increases, but does 
not become oscillatory or unstable. 

3. pitch angle takes care of itself within the system, 

| leaving the operator concerned only with the primary 
quantity, depth rate. 

4. the control system minimizes the effects of speed and 
hull geometry variations by increasing the natural 
frequency of the control loop by a factor of at least 
five. 

It is considered that this system, in conjunction with a pull-out 
time indicator discussed in the preceeding chapter, will provide near 
optimum depth changing and depth-keeping capability. Due to its 
broadband nature, one system of this type will provide adequate con- 
trol characteristics in a variety of hulls, obviating the need to 


tailor-design a control system for each new ship, 
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APPENDIX A 


COMPUTER STUDY 


A. Computer Simulation of the Proposed System 


The system as shown in Fig. IV was simulated on the analog 
computer in the following manner. To avoid loosing sight of the 
basic problem, the force and moment equations of the hydroframe as 
given in equations (3) and (4) were used to represent the hull. To 
facilitate the computer solution, these equations were rewritten in 


the following form 
d 


E O 
dr SE 


(3a) 100 Ge? 517 | 200 - 10 .1[2.' E (100 d0) + Kt (10 DE 


m' dr 


+ .1( 100 ' 外 











Sage _ Mala 
(48) 100 d%0 = .123 100w' = 2090 A „25 (1000 de), 
2 2 Cas dr 
dr U m 
+ .5 (100) + 45 Kt (10}(10 D 
where r= 10% 

2 
a = 10 = d = 100 É 5 
dz at” dr 


A point concerning the analog treatment of the lag network is 
worthy of note. Writing the defining relations of this component in 


the following form, 


Í 21 (s) E Ki (s + 2K) 
e (8 S + 2K 


and carrying through the indicated operations and change of time base 


yields 


ds de 
81 =Ki (10 주 ) + K (2K 4º = .2 


dr 


10 1) 
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Integrating once since = is unavailable, and noting that the 


initial conditions on the integration are zero yields 


10 844 = Ki (1010) + JIES e = 。2 91 | dr 
or Sy] 一 volK ] (100 e) + x | [-o2%,, (100 e) - ,23 dr 


Now the coefficient of the integral will appear in the computer 
diagram simply as a potentiometer for values of K between O and 10. 
K determines the location of the pole-zero pair of the lag network 
along the negative real axis of the s plane relative to the selected 
position at -2 and -.2, the ratio of their distances from the origin 
remaining constant. This facilitates a simple method of experimen- 
tally varying the phase shift introduced by the lag network without 
affecting the steady state error of the system. 


The feedback compensation is represented by the following: 


은 이겨 = ,310 (1000 afe ) + .163 (1000 de ) + .128 (100 8) - 


dr år” dr 


-.104 (10 E ) 


À limiter is employed to restrict the maximum value of 5, = dt 85 


to + 25º or .436 radian, 


dio is computed from the relation 
100 4" = 100 w' - 100 6 
go 


The complete computer diagram is shown on Plate I. 
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APPENDIX B 


ANALOG COMPUTER RESULTS 


Figs. VIII through XXV show the results of the analog computer 
study of the proposed system as plots of d 0, and e vs the non- 
dimensional time parameter # for various values of the three primary 

v 
parameters Ki» Z'/m', and U. The values of these parameters con- 
sidered are as follows: 
AS 
U Seo 520, 700, 9,0 


CO CEE A 


7 1 
W “1,0, 1,9, 2.0 


m' 


Curves for each value of U are shown on one plot for each pair of 


the parameters En and Ze considered, 





mi 
Fig. XXVI shows the performance of the system in pitch angle 


and depth rate for Zt /m! = -1.0, U= 5.0, K_' = 0.3, and K , = 9.25 


Figs. XXVII and XXVIII show this data for Ze of -1.5 and -2,0 


m! 


and all values of K considered with Ko adjusted to maintain the 


1 


product KA K, constant, 
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